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We study the equivalence of ensembles for stationary measures of interacting 
particle systems with two conserved quantities and unbounded local state space. 
The main motivation is a condensation transition in the zero-range process which 
has recently attracted attention. Establishing the equivalence of ensembles via 
• convergence in specific relative entropy, we derive the phase diagram for the con- 

densation transition, which can be understood in terms of the domain of grand- 
^ c*| canonical measures. Of particular interest, also from a mathematical point of 

Oh" view, are the convergence properties of the Gibbs free energy on the boundary of 

^ r~| that domain, involving large deviations and multivariate local limit theorems of 

subexponential distributions. 
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1 Introduction 

Zero-range processes are interacting particle systems with no restriction on the num- 
ber of particles per site, i.e. with unbounded local state space. The jump rate of each 
particle depends only on the number of particles at its departure site which leads to 
a simple product structure of the stationary measure U, Q]. These processes have 
recently attained much attention in the theoretical physics literature (see [3] and refer- 
ences therein) since they exhibit a condensation transition under certain conditions on 
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the jump rates |@]. If the particle density exceeds a critical value p c , the system phase 
separates into a homogeneous background with density p c and a condensate, where the 
excess particles accumulate. First rigorous results on a single species system [5] show 
that this phase transition can be understood mathematically in the context of the equiv- 
alence of ensembles. This is a classical problem of mathematical statistical mechanics 
|@] which arises naturally in the context of studying stationary measures of interacting 
particle systems with conserved quantities, such as energy or the number of particles. 
In general, interacting particle systems with several conservation laws are currently of 
particular interest in non-equilibrium statistical mechanics, since they show a very rich 
critical behaviour (see 10] and references therein). There are not many general results 
for such systems, for example for zero-range processes with more than one particle 
species there exist only non-rigorous case studies so far Jst]. The motivation of this 
paper is to understand the condensation transition in such multi-species processes on 
the rigorous level of the equivalence of ensembles. 

For simplicity of presentation we focus on systems with two conserved quantities, 
which we interpret as the number of particles in a two-species system. The local state 
space for each species is N = {0, 1, . . .}, i.e. the numbers of particles on each lattice 
site are unrestricted. We require that the process has a stationary measure of product 
form. Due to the conservation law this induces a family of stationary measures 7T£,n 
with fixed particle numbers N 6 N 2 on a lattice of size L, the canonical ensemble. 
Another standard family is the grand-canonical ensemble v^, where the numbers of 
particles are random variables. The densities p 6 (0, oo) 2 , the expected numbers 
of particles per site, are controlled by conjugate parameters p G M 2 , the chemical 
potentials. In our case vf^ is a product measure and is also defined for L — > oo, where 
we write z/ M . Let C K 2 denote the maximal domain, such that v x is normalizable 
with finite first moment. 

In the thermodynamic limit ~N L /L — > p as L — > oo with densities p, one expects 

that 



Bfi(p) e Dp : tt LiNl -> u^) as L -> oo . (1.1) 

The question of the equivalence of ensembles is for which values of p and in what 
sense (11.11) holds, and how it has to be modified in the presence of phase separation. 
The main results of this paper are: 

1. We establish the equivalence of ensembles (11.11) for all p e (0, oo) 2 under mild 
regularity assumptions on the stationary product measure. 

In the proof we use specific relative entropy (or relative information gain), which is 
based on results from information theory [9] and was previously applied to study large 
deviations and the equivalence of ensembles for Gibbsian random fields [10], marked 



point processes (|lll] and weakly dependent measures UM- A common feature of these 
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models is a bounded Hamiltonian, which corresponds to = IR 2 in the above setting. 
In this case, phase separation is a consequence of long-range correlations, leading 
to non-differentiability of the Gibbs free energy (or non-convex canonical entropy) 
and a first order transition fl3l In our case there are no spatial correlations, but 
typically C IR 2 due to the unbounded local state space, and condensation is a result 
of large deviation properties of v x on the boundary of D^, where it turns out to be 
subexponential. 

2. We show how the phase diagram for the condensation transition can be derived 
solely from the shape of D^, and explain its relation to the mode of convergence 
in (O) . 

The transition is continuous and is characterized by convergence properties of the 
Gibbs free energy on the boundary of D^. In the classification of [14] this corresponds 
to the case of partial equivalence of ensembles. 

Our results can be directly generalized to any number of particle species with ar- 
bitrary discrete local state spaces. We choose to work in a more specific setting for 
the simplicity of presentation, since it covers the basic novelties of the paper. From 
the point of view of non-equilibrium statistical mechanics these are the first rigorous 
results on the condensation transition in a system with several conservation laws. From 
a mathematical point of view, we adapt the theory of the equivalence of ensembles to 
study phase separation in systems with unbounded Hamiltonians. Even in the basic 
case of stationary product measures the different mathematical origin of the conden- 
sation transition leads to interesting new aspects. Our equivalence result involves a 
sharp condition on the number of particles and is valid on the (non-empty) boundary 



of D^. The analysis requires results on large deviation s 11151.11611 and multi-dimensional 
local limit theorems of subexponential distributions 1 17~ ll8|l . as well as convergence 



properties of multivariate power series similar to lfl9ll . In contrast to a previous study 
for single-species processes [5], the present paper provides a complete picture of the 
mechanism of condensation in a much more general context. 

Precise definitions of the ensembles and basic properties are given in the next sec- 
tion. The main results are given in Section 3, including the equivalence of ensembles 
and the construction of the phase diagram for the condensation transition. For com- 
pleteness, we also include some remarks on fluctuations and the spatial extension of 
the condensate (cf. [Hi|2^]). Proofs are given in Section 4. Since the main results apply 
for ensembles of measures in a general context, the paper up to this point is formulated 
without reference to zero-range processes, which are, however, the main motivation 
for this study. In Section 5 we explain why these processes provide a natural class 
of particle systems for the measures considered in the first sections, and illustrate the 
results on the phase diagram by several examples. Some results from convex analysis 
needed in the proof of the equivalence of ensembles are summarized in the appendix. 
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2 Preliminaries 



2.1 Canonical and grand-canonical measures 

Consider L independent identically distributed random vectors 

r){x) = (771 (x), 7/2(0:)) e N 2 , ieA t , (2.1) 

with some discrete index set A L of size |A L | = L. The state space X L = (N 2 ) Al 
is a measure space with <r-algebra induced by the product topology and the (a-priori) 
measure 

w L (r)) = J] w(rj(x)) G (0,oo) for 77 = {v(x)) xeAL • (2-2) 

xeA L 

This should be positive but not necessarily normalized, i.e. w : N 2 — > (0, 00) is 
arbitrary. Since X L is discrete, we simplify notation here and in the following by using 
the same symbols for a measure and its mass function, i.e. w L (rj) = w L ({ry}). 

We interpret the index set A L as a lattice of size L and 77 G X L as particle configu- 
rations of a two-species particle system. We do not specify the geometry of the lattice, 
boundary conditions or dynamics of this process, they should be such that w L is a sta- 
tionary weight, i.e. up to normalization, w L is a stationary distribution of the process. 
Generic particle systems with this property are zero-range processes discussed in Sec- 
tion 5. Apart from stationarity of w L , the only other requirement on the particle system 
is that the numbers of particles 

S L (ry) = (Si, E 2 ) (77) := £ r,(x) G N 2 (2.3) 

xeA L 

are conserved quantities for each species and that there are no other conservation laws. 
Then there exists a family of stationary probability measures h(H L ) w L which are ab- 
solutely continuous with respect to w L , where the Radon-Nikodym derivative depends 
only on the conserved quantities Ti L and can be written as a function h : N 2 — > [0, 00). 

The set of all stationary measures of the particle system is convex and the extremal 
measures are given by choosing h(H L ) oc 5s l ,n, i-e- proportional to the Kronecker 
delta, fixing the number of particles to N = (A r 1 , A^ 2 ) G N 2 . The family 

klAv) = J—U w(ri(x)) 5 Sl(t ,), n , NGN 2 (2.4) 

is the canonical ensemble and the measures concentrate on finite subsets 

X LiN = {r7|E L (77)=N} CX L (2.5) 



4 



of configurations with fixed particle numbers. The canonical partition function is 
Zl,n = w l (X L: n) E (0, oo), since 7r LjN = w L (. | {S L = N}) can be written as 
a conditional measure. By assumption, for each fixed L > 1 and N E N 2 the particle 
system is irreducible on Xlj* and ttl ; ts! is the unique stationary measure. All other 
stationary measures on Xl are convex combinations of canonical measures. 

Another generic choice is oc e M ' EL with parameters fi = (/ii,/U 2 ) E M 2 

called chemical potentials, defining the grand-canonical measures 

= -Az II w(fi{x))e^K (2.6) 

Each is supported on X^, i.e. is a random variable and the expected value is 
fixed by the chemical potentials fi, as is discussed below. These measures are particu- 
larly convenient since they are of product form. The normalizing (single site) partition 
function 

*(/*)= £>(k)e"- k (2.7) 

keN 2 

is an infinite sum, as opposed to models with bounded local state space, such as {0, 1} 
for lattice gases or { — 1, 1} for spin systems. For such systems, z(fj.) is defined for all 
fi E M 2 , whereas in our case the domain of definition of z will play a crucial role. 



2.2 Properties of grand-canonical measures 

We define 

L> M = {a* E R 2 J k i w ( k ) eM ' k < oo for i = 1, 2} . (2.8) 

keN 2 

This implies that for all /x E D^, z(fj,) < 00 and the product measure is well de- 
fined. Moreover, on the marginal v x has finite first moments, which are interpreted 
as particle densities and given by 

R = (Rx, R 2 ) : Dp -> (0, oo) 2 , where R^) = (77^1 , i = 1, 2 . (2.9) 

Here and in the following we write (..)„ for the expected value with respect to measure 
v. Note that Ri(fJt) = (77* (2?)) L independently of the lattice site x E Al, and that _D M 
as defined in (12.81) is the maximal domain of definition of R, i.e. = domR. We 
denote by 

Dp = R(.D M ) C (0,oo) 2 (2.10) 
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the range of R, which characterizes the set of all densities accessible by the grand- 
canonical ensemble. 

In the following we assume that w is exponentially bounded, i.e. 

3 £ £ (0, oo) Vk G N 2 : w(k) < £' k ' , (2.11) 

1/2 

where we write |k| = ||k||2 = (kf + fcf) ■ For convenience we further assume that 
the single site mass function w is actually defined on [0, oo) 2 with 

tofC 1 ([0, oo) 2 , (0, oo)) , (2.12) 

which imposes no restriction on the relevant values u>(k), k £ N 2 . 

Lemma 2.1 ^ (and thus D p ^ 0) if and only if < \2.11\) is fulfilled. In this case 
is convex and complete, i.e. 

A(/x*) := {/i | /ij < /i*, i = 1, 2} C D p whenever fi* £ . (2.13) 

Either = IR 2 or the boundary can be characterized in the rotated variables fi\ = 
/ii - /i 2 and jl 2 = 1^1+^2 by dD p = {(/}i,/2 2 (£ti))| Ai e ^}. Here ji 2 ■ R — > M w 
continuous and piecewise differentiate, with 

^2(^1) - - limsup (21ogw(k) + /^(A^ — k 2 )) / (h + fc 2 ) • (2.14) 

|k|^oo ^ ' 

A// ?/ze above properties also hold for domz, the maximal domain of definition of z. 
We have C dom z and int _D M = int dom z for the interior, so both sets are equal 
or differ only on the boundary. 

Since the grand-canonical measures are product measures the pressure is given by 
P(/*) = r lim t 1o ^ z ^ L = ' ( 2 - 15 ) 

L^oo Li 



which is the analogue of the Gibbs free energy. For all fi £ the density (12.91 ) can 
be written as 

Ri(fi) = d m p(») , i = l,2. (2.16) 

The derivatives are defined one-sided on dD^ n -D^,, which is possible due to complete- 
ness of Dp and the following lemma. 

Lemma 2.2 The single site marginal v 1 ^ has some finite exponential moments if and 
only if /x £ int D^. Moreover, p £ C°°(int D^, R), p £ C 1 (D Ai ,M) anJp « strictly 
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convex on D p . 

p and R can be extended continuously to 

fr- 00 Dp = { (-oo,^ 2 ) I G R : (jn,^) G D^} , (2.17) 

i.e. limits exist and are given by 

00 / \ 

p(-oo,/i 2 ) = ^w(0,A; 2 )e^ 2 and R(-oo,/i 2 )=( J . (2.18) 

k 2 =0 ^ °°i^)J 

An analogous result holds for d^'^D^. 

For i = 1,2, if D p is bounded in p it i.e. for all p G D p , pi < C for some C > 0, then 
D p is bounded in fa. 

Since p is strictly convex, R is invertible on due to (12.161) and we denote the 
inverse by M : D p — > D M . The entropy density s : (0, oo) 2 — > R of the grand- 
canonical measure (12.61) is the convex conjugate of the pressure given by the Legendre 
transform (cf. (IA.7I) ) 

s(p) = p*(p) = sup (p • p - p(/i)) . (2.19) 

Thus s, also known as the large deviation rate function, is strictly convex on D p and 
convex on (0, oo) 2 . For p G int-D p it is easy to see that p • p — p{p) has a local 
maximum at M(p) and thus 

s(p) = p • M(p) - p(M(p)) and M t (p) = <V(p) , i = 1, 2 . (2.20) 

Using convexity of D p and p(p) we can show that there exists a unique maximizer 
of the right hand side of (12.191) , also for p g" int _D p . This is the main result of this 
preliminary section. 

Proposition 2.3 For every p G (0, oo) 2 ?/zere exists a unique maximizer M(p) G -D^ 
of the right hand side of A2.19\) . such that 

s(p) = p-M(p)-p{M(p)). (2.21) 

M G C((0,oo) 2 ,R) and we have M(p) = M(p)for p G D p andM(p) G dD p f]D p 
for p £ D p . 

In particular, D p is closed in (0, 00) 2 an J <9-D p = R(<9.D p fl Z? M ), where dD p denotes 
the relative boundary of D p in (0, oo) 2 . 
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3 Main Results 



3.1 Equivalence of ensembles 

Consider a sequence of canonical measures ix L )Ni in the thermodynamic limit, i.e. 

~N L /L — > p as L — > oo with density p G (0, oo) 2 . (3.1) 

In the following we study the question if the sequence ^l,~n l converges to a grand- 
canonical product measure, and if yes, what is the mode of convergence. To quantify 
the distance between the measures we use the specific relative entropy 

^l,n(a0 = j^H(tt l ^\u^) , where 
^KN|^)=E^)log^M (3.2) 

is the usual relative entropy, since 7r i N is absolutely continuous with respect to v^. 
Using the relations 

vfa) z(p) L = w L (t]) e^ N for all rj G X L , N and 

vl ({S L = N}) z(p) L = Z LiN e^ N , (3.3) 

which are easily derived from (12.41) and (12.61) . we can write 



^,nO) = ~\ log^({S L = N}) =p(p) - ^ _ IlogZ L , N , (3.4) 

for all L > 1, N G N 2 and p G D^. 

The second part of (13.41) suggests that M(p) of Proposition l2.3l is the right chemical 
potential to minimize h,L,N L m the thermodynamic limit (13.11) . This is the content of 
the next theorem for which we need a further regularity assumption on the exponential 
tail of w, in addition to (12.1 II) and (12.121) . A convenient sufficient condition is that for 
all G [0, 7r/2] the limit in the radial direction 

lim - logw(reJ G R exists , (3.5) 

r— >oo r 

and is a continuous function of <\>. This can be relaxed considerably as is discussed after 
the proof in Section 4.3. (13.51) holds for example if w is convex, or if w = w\ + w 2 
where W\ is convex and w 2 has bounded derivative. 

Theorem 3.1 Assume A2.11\) . A2.12\) and ( 13.51) . Then for each particle density p G 
(0, oo) 2 and every sequence N L as in ( I3.il) 

lim h L ^ L (M(p)) =0. (3.6) 

L—>oo 
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From this result one can immediately deduce two standard formulations of the equiv- 
alence of ensembles, on the level of measures and on the level of thermodynamic 
functions. To formulate the first version we have to define all canonical and all grand- 
canonical measures on a common state space X = N A , where A is the (infinite) limit 
lattice of an appropriate sequence (Al)z=i,2,..- The precise construction is deferred to 
Appendix B, since it is only necessary to formulate (13.71) and has no further importance 
for our results. 



Corollary 3.2 For each p G (0, oo) 2 we have 

{f)*L,n L -» (f)»u ip) asL^oo, (3.7) 

for all cylinder test functions f G C(X,M.) with {^)u M ,- ) < oo for some e > 0. In 
particular, this includes all bounded f G Cf,(X, M), which is equivalent to convergence 
in distribution. Moreover, 

lim — logZ L>Ni = -s(p) . (3.8) 

For p G int D p , has some finite exponential moments by Lemma |2~TI so in par- 
ticular the corollary implies convergence of the local densities f(rj) = r)i(x). We note 
that for a single species with p G int D p convergence is shown even for L? test func- 



tions in plQ . Appendix 2.1. The proof given there relies on rather involved estimates 
on the rate of convergence in the local limit theorem, whereas the proof via relative 
entropy is much simpler (see section 4.3). Moreover, our result covers several particle 
species and can be generalized to p G" D p , which is the main point of this paper. In this 
case the nature of the convergence changes and (13.71) is violated for f(rj) = r)i(x) for at 
least one species i, as will become clear in the next subsection. This difference in the 
mode of convergence is a result of the unbounded local state space and is a signature 
of the condensation transition. 

For systems with bounded local state space (13.61) implies convergence for all cylin- 
der test functions / G C (X, R) . But in case of p G" D p the limiting measure would be 
a mixture of grand-canonical measures, corresponding to coexisting domains with dif- 



ferent distributions for large finite systems (see e.g. 1130). This phenomenon is called 
phase separation. In analogy to this classical case we interpret our limit result in the 
following way: For p G" D p the system phase separates into a (homogeneous) back- 
ground phase with product measure t^( p -\ given by Theorem 13 . 1 1 and a condensate 
or condensed phase which contains the excess particles. According to (|3.7I) . the con- 
densate cannot be tested by cylinder functions in the infinite system, its existence is 
only a consequence of the conservation law (in contrast to classical phase separation). 
The interpretation for large finite systems is that the volume fraction covered by the 
condensate domain vanishes as L — > oo. In fact, this domain typically concentrates 



9 



only on a single lattice site, which is proved under additional assumptions in Section 
3.3. 

Due to the conservation laws, the phase space of the particle system is (0, oo) 2 , the 
set of densities p. We say that the particle system exhibits a condensation transition, if 
D p C (0, oo) 2 . As order parameter of the phase transition we choose the mapping 

R c : (0, oo) 2 -> D p , with R c (p) := R(M(p)) . (3.9) 

According to the above interpretation, R c (p) describes the density of the background 
phase in a system with global density p. Note that by Proposition [23] 

so H c (D p ) = D p and R c is a projection from (0, oo) 2 onto D p . By Lemma l2~2l and 
Proposition 12.31 R c G C((0, oo) 2 , Dp) so the transition is continuous (second or- 
der), which is directly related to the fact that p G C 1 (Z^ M ,1R). This is in contrast to 
systems with bounded local state space, where we would have Dp = R 2 and non- 
differentiability of p would lead to a first order phase transition with discontinuous 
order parameter [[lj, [HI] . 

3.2 Phase diagram 

In this section we apply standard results from convex analysis, which are summarized 
in Appendix A, to characterize the phase diagram of the system. By Proposition I2.3L 
dDp = H.(dDp n Dp), and thus condensation occurs if and only if dDp n Dp ^ 0. 

Theorem 3.3 For every p c G dD p with p = M(p c ), the preimage R~ 1 (p c ) is given 
by the subgradient Sp(/j,) as defined in < \A.6\) . Moreover, 

/a / {p c + An/,) A > 0} , dDpdiff'ablein p 

OPW ~ \{ Pc+ A+n++A-n- | A + ,A" > 0} , otherwise ' ^' iiJ 

n M denotes the normal vector to dDp in p and n^, n~ the two limiting normal vectors, 
in case dDp is not differentiable in p. 

Note that by convexity of Dp, n+ and n~ are well defined as the extremal normal 
directions to the set of supporting hyperplanes in p. In case that the points of non- 
differentiability of dDp (see Lemma [27Tb accumulate in p, 11+ = is also possible. 
In the following we use Theorem l3.3l to construct the phase diagram and its properties. 



By definition (13.91) , the preimage 

Rc\Pc) = {P e (0, oo) 2 | R(p) = p c } (3.12) 
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denotes the set of all densities having background density p c E dD p , and (13.1 II) implies 
that this is a linear set in normal direction to D M . By Lemma |2~T1 the normal vectors to 
dD p have two nonnegative components. So a first direct consequence of Theorem [33] 
is that (as expected) 

R c (p)<P i.e. R c ,i{p)<Pi for« = 1,2. (3.13) 
We say that species i condenses if R c ,i(p) < Pi an d define 

A i = {pe (0, oo) 2 I R Cji {p) < Pi }c (0, oo) 2 \ D p . (3.14) 
So the phase space (0, oo) 2 can be partitioned in the following way: 

• homogeneous phase region D p , R c (p) = p, 

• condensed phase region A\ \A 2 , condensation of species 1 only, 

• condensed phase region A 2 \ A\, condensation of species 2 only, 

• condensed phase region A\ fl A 2 condensation of both species, 
defining the phase diagram of the model, 

PD = [D p , A 1 \ A 2 , A 2 \A 1 ,A 1 nA 2 }. (3.15) 

The topology of the phase regions depends on the weight w and all phases except D p 
may also be empty. Examples are given below in Figures Q] and [2] and in more detail 
also in Section 5. 

The entropy density s and the pressure p are convex conjugates (12.191) , and thus 
p E 5p(p) -v=> p E ^ s (p) ( see Theorem I A.3I) . Together with Theorem [33] this implies 
that s is an affine function on R c 7 1 (p c ) for every p c E dD p , i.e. 

s(p) = s(p c ) + (p-p c )M(p c ) forallpGR; 1 ^). (3.16) 

So s has a non-strictly supporting hyperplane in p c , i.e. we have partial equivalence of 
ensembles in the sense of I bill . 

By M every condensed phase region is mapped on dD^ fl D^ L . Due to their def- 
inition, A\ \ A 2 and A 2 \ A\ correspond to special parts of that boundary. Define 



d i D il :={pE dDp \ dDp diff'able in p, n M || ej for i = 1, 2 , (3.17) 

where ej denotes the unit vector in direction i. By regularity of Dp given in Lemma 
12.11 namely convexity and completeness, d l Dp (if non-empty) is a straight line of the 
form 

d 1 Dp = [p 1 + Ae 2 | A E R} or d 1 D fl = [p 1 - Ae 2 1 A > 0} , (3.18) 
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for some /x 1 G dD^. In the first case d 2 D ll has to be empty, and an analogous version 
of (13.181) holds for d 2 D ll . The connection with the phase region is basically a direct 
consequence of Theorem 13 .31 summarized in the following. 

Corollary 3.4 The phase regions D p , A\ \ A 2 and A 2 \ A\ are simply connected, 
M(Ai \ Aj) = d { D p n Dfj, and 

A i \A 3 =R(d i D fM nD fl ) + {Xe i \\>0} for i = 1, 2 , j ^ i . (3.19) 

Here we use the convention that M\ + M 2 = {mi + m 2 \rrii G Mi} which is empty if 
one of the sets is empty. 

D p 7^ and it is connected to (0, 0). If Ai \ Aj ^ 0, it is connected to {p \ pj = 0}. 

Using this, the phase diagram is uniquely specified by dD^ up to the exact location 
of dD p = R(<9/} M fl Dfj). The phase region A\ fl A 2 corresponds to a more complicated 
part of the boundary. It may in general be disconnected and is most easily found as the 
complement of the union of all other phase regions. As a direct consequence of (13.191) . 
the phase boundary between A\ fl A 2 and A { \ Aj, if both are non-empty, is a straight 
line in direction ej. These results are illustrated in Figured] and [2] for 

^( k ) = Tx , i.\ i — r~o ' ( 3 - 2 °) 



[l + b) kl \h + 2 

where + = n^=o 1 (*"'"^) denotes the Pochhammer symbol. This is the stationary 



weight of a zero-range process that has been introduced in [|8|], where also the phase 
diagrams are derived. This derivation and the corresponding zero-range process will 
be revisited in Section 5. 

By Proposition 12.31 occurrence of condensation can be characterized by boundary 
properties of D^. For single species systems this can be directly translated to a condi- 
tion on the weight w. dD^ = {fi c } is only a single point and /i c G if and only if 

i/J o (Jfe) oc w{k) e^ k = o{k- 2 ) , (3.21) 

and v x has finite first moment p c , the critical density. For a two species system the 
condition dD^DD^ ^ cannot be rephrased as a simple condition on w, and in general 
knowledge of is required. But for specific systems this is usually not difficult to 
obtain, and combining the results of this section, this is also sufficient to derive the 
phase diagram. This is explained in more detail for some generic examples in Section 
5. 

3.3 Properties of the condensed phase 

Note that the canonical measure can be written in conditional form 

ttl,n £ =it(.\{Z L = N L }) (3.22) 
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Figure 1: and phase diagram for w as given in (13.201 ) with 6 = 4. Dashed lines on the left 
denote the normal directions to dD^ and on the right they determine the function R c (p) and 
the phase regions as given in Corollary 13.41 



for all p G D^, and in particular for p = M(p) where p = lim N L /L. For p £ D p , 

L^oo 

{H L = N L } is a large deviation event for v^,y A typical configuration for 7r i Ni 
then corresponds to a most likely event to realize this large deviation. 

Theorem 13.11 implies that in such a configuration for large L we have coexis- 
tence of a critical background domain with density R c (p) and a condensate containing 
(p — R c (p))L particles on average. The number of particles in the condensate fluctu- 
ates, and due to the conservation law the distribution is given by the fluctuations of the 
number of particles in the background domain. The latter has extensive volume and the 
distribution converges to a product measure, so the particle number is approximately 
given by a sum of iidrv's. Thus if has finite second moment the central limit the- 
orem applies and the fluctuations should be Gaussian. If the second moment is infinite 



the fluctuations are non-Gaussian and determined by the Levy stable law of u^r ^ [1 1 VM . 
This picture is consistent with results on single- species zero-range processes Q22Q. In 
general, the background distribution for each condensing species is subexponential, so 
non-Gaussian fluctuations are indeed possible. 

Lemma 3.5 Suppose Ai ^ 0. The single site marginal v]l {k) = ^2 ^(k) with 

It* Ayj 

respect to species i has a subexponential tail, i.e. 

liminf-^logz/^fA;) = , (3.23) 

k~>oo k 

if and only if p £ M(Aj) C dD^ H D^. This holds with lim instead of 'lim inf ifw has 
a regular tail as given in di.51) . In this case there also exists a subexponential sequence 
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Figure 2: and phase diagram for w as given in (13.201) with b < 3, which implies A\ n ^2 = 
(for the relevant calculations see Section 5). For 6 = 3 phase A\ \ A2 is non-empty as shown 
on the left, for 6 = 2 only species 2 condenses as shown on the right. 

in directions normal to dD^. More precisely, for all n M normal to dD^ at /x and for 
every sequence k n with |k n | ^00 and k n /|k n | — > n M as n — > 00 we have 

lim --Llog^(k n ) = 0. (3.24) 
J k n I 

The statement (13.241) is particularly important for the proof of Theorem 13.11 It is re- 
markable that a subexponential sequence exists precisely in the direction normal to 
dD^. Together with (13.111) this allows all the excess mass in the system to condense 
on a single lattice site. If dD^ is not differentiable in fi, (13.241) should also hold in 
all intermediate directions between n+ and n^. This is supported by the examples we 
studied so far (some of which are given in Section 5), which all show condensation 
on a single lattice site. However we are not able to prove this in general and it is also 
not required to prove Theorem 13. 1[ Statements similar to (13.241) on the limit behaviour 
of indices of multivariate generating functions have been derived only recently in the 
combinatorics literature based on Cauchy's integral formula (see dl9D and references 
therein). 

'if dD^ is differentiable in fi, is unique. Otherwise the statement holds for the two limiting 
normal directions n+ and (cf. Theorem l3.3b . 
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Regarding the spatial extension of the condensate domain in large finite systems of 
size L, Theorem 13.31 only assures that it covers a non-extensive volume of o(L) sites. 
Suppose the tail of the single site marginal with respect to a condensing species is a 
power law with finite first moment, i.e. 

(*) ~ k- b for some b > 2 as k -> oo . (3.25) 

Then the following statement holds. 

Theorem 3.6 For p e A{, p = M(p) and with t\3.25\l we have a weak law of large 
numbers, i.e. for all e > 



maxr]i(x) - - R c ,i(p)) 



> e) ^ as L — > oo , (3.26) 



w/zere z>£' W = z/£' W ( . | {S' L = iV^}) w*/z iV JiL /L -> Pi . 

For single- species systems vr ii Ar L = z/^(. |{£ L = A^}) and the theorem implies 
convergence with respect to the canonical distribution til.n l iSl- Therefore in typical 
configurations for large L all the excess mass concentrates on the site with maximal oc- 
cupation number, so the condensate covers only a single lattice site. Since the particle 
system is translation invariant, the location of this site is chosen uniformly at random. 
This property i s ty pical for large deviations of distributions with subexponential tails. 
The results in 11151 llol and Monte-Carlo simulations of zero-range processes suggest 
that Theorem 13 .61 also holds for other subexponential tails. The proof given in Section 
4.4, however, which is a slight extension of a result in Il20ll . works only for power law 
tails. In the case of two or more species, til,n l 

7>4' (i) (- I {^l = N i,L})- However, 
Monte-Carlo simulations of zero-range processes |23j] confirm that the above implica- 
tion of Theorem l3~6l which is rigorous for a single species, also holds for two-species 
systems. 



4 Proofs 

4.1 Proofs of Section 2 

Proof of Lemma I2TT1 Note that = domR = domi?! fl domi? 2 and 

By definition, z(p) > so domi?j is the domain of convergence of the sum in the 
numerator. Suppose domi?j is non-empty and contains p 1 ^ p 2 . For each p = 
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A/it 1 + (1 — A)/x 2 , A E (0, 1), by convexity of the exponential, 

e M-k < A g^ k + (i - A) e^ k (4.2) 
for all k G N 2 . So each term in the sum (14.11) can be bounded and we get 

kiw{k)e* w < A J2 ^(k)e Ml k + (1-A) ^ hw^e^* < oo , (4.3) 

kgN 2 keN 2 keN 2 

so fj, E domi?j. Thus is the intersection of two convex sets and is itself convex. 
Since ^] keN 2 h w(k) e M ' k is monotonically increasing in ii x and fi 2 , /^* £ implies 
H E for all // E A(/x*) and is complete. Replacing /cjU'(k) by w(k), both 
arguments also hold for domz, the domain of convergence of XlkeN 2 W S -) e^. So 
dom z is convex and complete, and certainly C dom 2;. 

To get a representation of the boundary we change variables and write 

00 

J R.( /X ) = J2 Yl ki w ^> e^-^^e^ 2 , (4.4) 

m=0k 1 +k 2 =m 

with /i! = /i! — yu 2 and /i 2 = /ii + /i 2 . The boundary of .D^ is then given by 

ft 2 {fi 1 )/2 = -limsup — log( ^M^e^ 1 ^ 2 ) . (4.5) 

Since for every positive function / 

sup /(k)< ^ /(k)<(m + l) sup /(k), (4.6) 

ki+k2=m , , ki+k2=m 

ki+k2=m 

and the logarithm is monotone we have 

/^(Ai) = _ limsup ( 21ogw(k) + /^(fci - fc 2 ) )/ (&i + fc 2 ) . (4.7) 

|k|^oo V ' 

This is independent of i = 1,2 because logki/(ki + k 2 ) — ► 0, proving (12.141) . In 
particular, replacing &^(k) by w(k) leads to the same formula for dom z. So and 
dom z can only differ on the boundary, which implies int dom z = int D^. 

Since is convex, jl 2 is piecewise differentiable (see e.g. [24], Chapter 1). By 
completeness of it is clear that 7^ if and only if /i 2 (0) > —00, which is 
equivalent to 

limsup — logw(k) = inf sup — logw(k) < 00 . (4.8) 

|k|^oo h + k 2 neN | k |>„ ki + k 2 
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This in turn is equivalent to (12.1 II) . because if (12.1 II) holds (14.81) is bounded by £ < oo, 
and if (|2.1 II) does not hold there exists a sequence k n for which log w(k n ) is un- 
bounded. The same argument works directly for dom z. □ 

Proof of Lemma 12.21 The existence of some finite exponential moments follows 
from the identity 

(e 9 -*>}„ h = Z(n + 0)/Z(n) (4.9) 

which is finite for sufficiently small E M 2 if and only if /x G int D^. Therefore 
p E C°°(int D^, E) and the covariance matrix 

D 2 p(») = DBJja) = Cov(^) = «M>y iJ=1>2 , (4.10) 

is well defined on int D^, where (rji r\j) c vl := (rji T)j) v i — (r)i) v i (Vj)^ • Further, D 2 p(^i) 
is symmetric and positive definite, because 

a T -(D 2 p(M))a=<(a-(r /l ,r /2 )) 2 )^ 1 >0 (4.11) 

for all a E M 2 with |a| = 1. Hence the eigenvalues of D 2 p(n) are real and positive, 
which ensures that p is strictly convex and R is invertible on int D^. 

Moreover, by completeness of D^, one-sided derivatives d^p exist on dD^ D D^, 
so p E C 1 (£) /i , E). Strict convexity of p extends to dD^ fl D^, since for any linear 
subset of the boundary + Ae^, in direction e^, parametrized by A, 

d\p(n + Ae^) = ■ R(/x + Ae^,) is strictly increasing with A . (4.12) 

This holds because by completeness of D^, E [n/2, 7r] (or equivalently [37r/2, 2n]), 
and for E (n/2, 7r] and G [7r/2, tt), — i?i and -R 2 , respectively, are strictly increasing 
with A along + Ae^. 

With a change of variables ipi = e Mi G (0, oo), i — 1, 2 we get 

z(/i) = 5(0) = ^ «;(k) ^ ^ • (4-13) 

keN 2 

The domain of convergence of this two dimensional power series certainly contains 
points with ipi = or ip 2 = 0, corresponding to [L\ = — oo or // 2 = — oo, respectively. 
So the definition of z can be extended to d % ~ cx 'D ll , i — 1, 2, and thus also the definition 
of p, as 1(0) > W(0, 0) > for all 4>. Since z is continuous on its domain of 
convergence and fel = 5 kl:0 , 

oo 

lim z(0) = z(0, ip 2 ) = V io(0, fc 2 ) </> 2 2 (4-14) 

0^(-oo,^ 2 ) ^— ' 

^2=0 
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which yields the result for p(—oo,fi 2 ), an d analogously p(fi\, — oo). An analogous 
argument holds for R, where 



- oo 

i? 2 (-oo,/i 2 ) = Yh w (0,k 2 )e^ k2 (4.15) 

z[— oo, u 2 ^— ' 

can be written as i? 2 (— oo, /i 2 ) = dp 2 p(—oo, fi 2 ), and oo, yU 2 ) = for all /i 2 , since 
^0 fcl = for all h. 

By completeness of Dp, there exists /i 2 G (— oo, oo], such that d l ~°°Dp = {— oo}x 
[— oo, or {—00} x [—00, fji^], where the second case is only possible if ^ < 00. 
Again by completeness, Dp is bounded in /i 2 if and only if p* 2 < 00. 
Now we can proceed analogous to ||2lll . Chapter 2, Lemma 3.3: Let Dp be unbounded 

in fi 2 , i.e. d l ~°°Dp = {-00} x [-00, 00). Suppose R 2 (-oo, fi 2 ) = dp 2 p(-oo, fi 2 ) < 
C for all fi 2 G [—00, 00). Then 

log Z{ ~°°' ^ = p(-oo, fi 2 ) - p(-oo, 0) < C// 2 , (4.16) 
z(— 00, (J) 

and thus z(— 00, fx 2 ) < z(— 00, 0) e CM2 as fi 2 — > 00. This is in contradiction to 2: being 
a power series with positive coefficients, and so R 2 is unbounded. Since R 2 is also 
monotone increasing on d 1 ~°°Dp this implies lim R 2 (—oo, /x 2 ) = 00, and D p is un- 

bounded in p 2 . □ 

Proof of Proposition l2.31 The proof uses basic results of convex analysis which are 
summarized in Appendix A. For fixed p G (0, oo) 2 , the function F : R 2 — ► (—00, 00], 

F(, ^JfW-P'M,^dom, 
v 7 [ 00 , p ^ domz 

is strictly convex on Dp C domz. So by Theorem IA.ll argmhiF R the set of 
minimizers on Dp contains at most one element. By (12.191) . s(p) = — inf Mg £> M F(p), 
so M(p) is a minimizer of F. Also F G C 1 (D A1 , R) (by regularity of p), and 

VF(p) = p - R(p) . (4.18) 

Thus for p<ED p we have VF(M(p)) = and F has a local minimum in M(p). By 
Theorem lA.il this is then the unique global minimum of F on Dp. On the other hand, 
if p D p , VF(p) 7^ (0, 0) for all p G Dp, and if there exists a minimizer for F it 
has to be in <9dom z. In the following we will show that this minimizer exists and lies 
in fact in dDp n Dp. 
First we show that 

inf F(p) = inf F(p) for some compact, non-empty E C Dp . (4.19) 
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Since F is convex and C 1 we have for a fixed p* G int _D M 

F(p) > F(p*) + VF( M *) ■ (p - A**) = + " p) ' (/* " A**) , (4-20) 

for all /x G M 2 , following (IA.6I) . For any p G (0, oo) 2 , there exists p* G int .D^ such 
that Ri(fj,*) < pi, i = 1,2, since R(p) — >■ (0,0) as p — > (— oo, — oo) and is 
complete and nonempty. Thus 

F(n) > F(p*) for all p G -E* := {p : (R(p*) - p) ■ (/x - /x*) > 0} . (4.21) 

Since both components of R(/x*) — p are negative, the half-space E* contains the cone 
A(/x*), as given in (12.131) . The boundary dE* is a straight line. Since p* G int domz 
and dom z is convex, dE* fl <9dom 2 consists of at most two points. 

1. Suppose dE* n <9dom z has two points, then dom z \ E* is bounded and contains 
the level set E :— \ev< F i^F. By definition (14.171) . F is lower semicontinuous 
and thus E is closed by Theorem IA.2I and therefore compact, which implies 
(EH. 

2. dE* fl <9dom,2 = is only possible if ddomz || <9-E* is itself a line. By strict 
convexity of F the level line lev =J F( M *)F touches c^* only in the single point 
p*. So we can find a point ji** close to p* such that dE** fl <9dom z contains one 
point, where E** is defined as in (|4.21l) . Then continue as in case (3) below. 

3. If dE* fl <9dom z consists of one point, dE* fl dom z is a semi-infinite line and 
let be the direction of that line in which there is no intersection point. Then 
<p G (7r/2, 7r) or G (37r/2, 27r), and in the first case dom 2 is unbounded in p 2 
and in the second case it is unbounded in Concentrating on the second case, 
by Lemma l2~2l there exists p G such that i?i(p) > p\, and thus by continuity 
of F and (14.181) there has to be a point p 1 G such that d fll F(fi 1 ) > and 
d fia F(fjt 1 ) < 0. Then domz \ (E 1 * U E 1 ) is bounded, where is defined 
analogously to E* in (14.211) . If p 1 G E* then £ := lev< F(/J .)F C domz \ 
(£■* LIE 1 ) is bounded and compact as in the first case (1). If p 1 ^ E* we 
take E :— lev< F ( M i)F C dom z \ (E* U E 1 ), which is also compact. The 
same argument works for <\> G (ir/2, it), with p 2 such that d fll F(fi 2 ) < and 
^ 2 F(p 2 ) > 0. 

This completes the proof of (14. 1 9b . 

Since F is continuous on dom z, (14.191 ) implies that it has at least one minimizer 
p on domz. By completeness of domz and D^, p + te G for all t < 0, where 
e = (1, l)/y/2. By strict convexity of p, 

— p(p + te) = R(p + te) • e (4.22) 
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is monotone increasing with t. Therefore lim t ^ R(p + te) ■ e exists and is given by 
lim ]T (fcx + fca) e^ +te )' k = R 1 (fi) + R 2 {p) (4.23) 

by monotone convergence. Now suppose p £ Dp, then this implies 

d F(p + te) = (R(p + te) - p) ■ e -> oo as t / . (4.24) 



This is a contradiction to p being a minimizer of F, so /i G Dp. By the first part of the 
proof we have uniqueness on D p , so M(p) := p G <9Dp fl D p is the unique minimizer 
of F for all p £ D p . 

Since F G C l (D p , M) is linear in p and strictly convex on D p for all p G (0, oo) 2 , 
its maximizer M(p) is a continuous function of p. Since R and its inverse are contin- 
uous, D p and D p are diffeomorphic and dD p D D p = H(dD p fl Dp). Suppose there 
exists p G <9D p \ -Dp. Then there exists pj G D p ,j G N, such that p — > p as j — ► oo. 
But then by continuity of M and R 

Pi = R(M( Pj )) - R(M(p)) G D p , (4.25) 
in contradiction to p ^ D p . Thus D p is closed in (0, oo) 2 and <9D p = R(<9Dp fl D M ).D 



4.2 Proofs of Section 3.2 

Proof of Theorem 13.31 In the following suppose that dD p is differentiable in p c . In 
case it is not, analogous arguments hold with An replaced by A + n + + A~n~. 

For every p c G dD p also p c G D p since D p is closed, and p c = M(p c ) G dD p nD p 
by Proposition l2.3[ One-sided derivatives of p exist in p c and 

Vp(p c ) = p c G 5p(p c ) , (4.26) 

where the definition of the subgradient is given in (IA.6I) . Let n _L dD p in p c . Then 
An ■ (p — p c ) < for all A > and p G Dp, and thus 

p(p) > p(p c ) + (p c + An) • (p - p c ) , (4.27) 

i.e. (IA.6I) holds and p c + An G 5p(p c ). On the other hand, if n / <9D P or A < there 
exists p* G D p with An • (p* — p c ) > and p(p*) > p(p c ) (by strict convexity of p) 
and (IA.6I) does not hold, so that 5p(p c ) = {p c + An | A > 0}. 
p c = M(p c ) minimizes p(p) — p c ■ p and thus minimizes also 

Km) - Pc ■ P ~ ^ n ■ (A* - A* c ) > K/0 ~ Pc ' ( 4 -28) 
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for all A > 0, with equality if and only if p = p c . So M(<5p(p c )) = {p c }, 
K c {5p(iJb c )) = {p c } and thus 5p(p c ) C R~ x (p c ). On the other hand, p G R c 7 1 (p c ) = 

M (p c ) implies M(p) = p c , and thus s(p) = p ■ p c — p(p c ) by Proposition I2.3L 
which is equivalent to p G Sp(p c ) by Theorem IA.3I Therefore 8p(p c ) = R~ 1 (p c ), 
finishing the proof. □ 

Proof of Corollary 13.41 Since {p|pi < p 2 (0)/2, i = 1,2} c by the proof of 
Lemma |2~T1 and R(p, //) — >■ (0, 0) as p — >■ — oo, D p ^ is connected to (0, 0). D p is 
simply connected since it is diffeomorphic to the convex set D p . 

By definition of Aj \ Aj in (13.141) and Theorem 13 .31 for each p G A; \ Aj, e, is 
a normal vector to dD p at M(p), and thus M(p) G n D M by definition (13.171) . 

On the other hand, for every p G d*/^ n D p , M _1 (p) D R(p) + {Ae. ; |A > 0} by 
Theorem 13 .31 and thus we have shown (13.191) . Since d l D p is either empty or a (simply 
connected) straight line according to (13.181) . the same holds for d l D p n D p by com- 
pleteness of and thus also Aj \ Aj is simply connected. If d l D p R -D M is non-empty 
it is connected to \ij = — oo by (13.181) for j ^ i, and thus H(d l D p fl D p ) is connected 
to {p\pj = 0} since Rj(p) — > as p.,- — > — oo. Therefore by (13.191 ) shown above, 
\ Aj is connected to {p\p ] ■ = 0}. □ 



4.3 Proof the equivalence of ensembles 

In the proof of Theorem l3.1l we make use of Lemma [331 which is proved first. 

Proof of Lemma 1331 If a = lim inf — — log v 1 '^ (k) > then for almost every k 

k^oo k 

it is v]l {€) {k) < e~ ak . Thus 

, ,„. z(u + eia/2) 

(e* a /2) = ^ . ; ' ' < oo , (4.29) 
N 1 v » z{fjt) 

which is in contradiction to p G M(AA according to Lemma |2~T1 

To prove the second statement of the lemma we use the integral criterion that for 
all p G M 2 , z(p) < oo if and only if 

p poo /*7r/2 

z(p)= w(*)e^d 2 x = j I w(re^)e r ^rd<pdr < oo . (4.30) 
J[o,od) 2 Jo Jo 

By regularity of w we can write the integral in polar coordinates (r, <p) with r = ||x|| 2 
and = x/r. Depending on w, (14.301) is not necessarily equivalent to < oo. 
In this case, however, w can be replaced by another function of the same regularity 
having the same values w(k), k G N 2 , for which equivalence holds. For example one 
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could take w to be a smoothed version of a piecewise linear interpolation of the points 

a;(k),kG N 2 . 

First assume that dD^ is differentiable at /x G dD^ and pick a direction G 
[0, 7r/2]. Define the positive half space 

Q(/x,0) = {6>|(0-/x)-e >O}cM 2 . (4.31) 

Since Z)^ is convex and dD^ is differentiable in /x we have 

Q(/x, 0) n int = if and only if = n M 1 <9£> M . (4.32) 

Let ^ €= [0, 7r/2] be the direction of n M . Then for all ^ -0 there exists /x^ G Q(/x, 0) n 
int .D M , and 

w(re^) e r/ ^ = w(re^) e r/ V e ^ e^^* . (4.33) 

Since i/^ has some exponential moments and (fx — (jt^) ■ < 0, there exists > 
such that 

w(re^) e r ^ = 0(e~^ r ) . (4.34) 

So the integrand of (14.301) decays exponentially fast in all directions ^ ip. Further, 
for all G [0, 7r/2] the maximal to choose is given by 

Za, = - lim - logw(re^) - e^ - /x , (4.35) 

r— »oo f 

The limit exists due to regularity of w and further is a continuous function of 0. We 
know that > for all 7^ ip an d thus in direction ip the integrand of (14.301) can 
decay only subexponentially, i.e. = 0. Otherwise there would exist with 6>.j > 
such that 6 • < e^/2 for all G [0, 7r/2]. By regularity of the integrand we could 
use Fubini's theorem to get 

n-K/2 1 

z(n + 0)< / d0<oo, (4.36) 

Jo e 4> ~ " ' e 4> 

which is in contradiction to /x G dD^. Alternatively, we also have — > as — >■ -0, 
since 

sup |0 — /x| — >• as — ► -0 , (4.37) 

9£Q( f j,,<p)nD ti 

and e^, = follows by continuity of e^. 

This second argument also works if dD^ is not differentiable in /x for the two lim- 
iting normal directions ip~ < ^P + - (14.371) holds as /* ip~ and \ ip + , leading to 
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e^p- = and = 0, respectively. With continuity of w the statement for u l ^{k. n ) 
follows for every sequence k n with |k„| — > oo and k„/|k n | — > n M . □ 

Proof of Theorem 13.11 According to (13.41) the specific relative entropy is 

h L ^ L (M(p)) = -1 log4 (p) ({S L = N L }) , (4.38) 

and we have to find a subexponential lower bound to the probability on the right-hand 
side. 

For p G int D p we have M(p) = M(p) G and ^m(p) nas exponential moments 
and thus finite covariance. The limit distribution of (£^(77) — pL) j\[L is a bivariate 
Gaussian and we have ^/ p )({^i(*7) = N^}) ~l/LasL— >ooby the multivariate 
local limit theorem, since Nl/L — > p. For p G dD p the same argument holds as 
long as ^M(p) nas fi^te covariance. In case it does not, the same conclusion follows 
from convergence to an a-stable limit law with rate of convergence 1 / L 2 / a f(L), where 
a G (0, 2) and / is slowly varying as L — > 00 (see [Q/zLQil] for multivariate local limit 
theorems in the non-Gaussian case). 

For p G" D p a typical number of particles under v^r s is [Rc(p)£], where [..] 
denotes the integer value. We get a lower bound for ^537 \ ({^l = Nl}) by a special 
configuration where the excess particles are concentrated on the first lattice site, i.e. 
77(1) = N L - [R c (p)L] . Thus we have for all p G D p 

^({£ L = N L }) >^({S L = N L ,77(1) = N L - [R c (p)L]}) . (4.39) 

First we assume that <9.D M is differentiable in M(p). By Theorem 13 .31 p — R c (p) is 
perpendicular to dD p in M(p). Thus by Lemma [331 

fc L|NL (M(p)) < - I log (N L - [R c (p)L] ) 

-Z l0g ^)({ S ^ = W^]}) ^ ( 4 - 4 °) 

as L — > 00, since (N^ — [R c (p)L]) /L —> (p — R c (p)). The second term vanishes 
with a local limit theorem analogous to the case p G <9-D p . If <9£> M is not differentiable 
in M(p) we may have two limiting normal vectors n ^ p ) an d n ^j( p ) which are linearly 
independent. Therefore there exist a\, a~[ > such that 

N L - [R c (p)L] = [a+n+ J + [e^nZ_J (4.41) 

for all L. Then we proceed analogous to (14.401) and get a lower bound by distributing 
the excess particles on the first two lattice sites according to 77(1) = [^nt ^] and 
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ry(2) = [ a L n M( p )] • By Lemma 1331 we have as L — ► oo, 

h L ^ L (M(p)) < - i log *4 (p) ( [a+n+ J ) - - log ( [o^mJ ) 

-I lo s^p)({ s ^ = [ R ^) L ] }) - • < 4 - 42 ) 

□ 

We note that condition (13.51) is not necessary to conclude (14.401) and (14.421) . It would 
be enough to have some sequence k n in direction p — R c (p) such that (13.241) holds in 
the limit n — > oo and sup n k n+ i — k„| < C for some Cel. However this is still 
not a necessary condition and therefore we used the comparatively simple assumption 
(13.51) which is fulfilled by a large number of examples. 

Proof of Corollary HH The first statement (Q with / e C(X, R) such that 
(e^)^ < oo for some e > is shown in ||2~5"ll . Lemma 3.1. It follows directly from 

subadditivity of relative entropy (see e.g. [§]) and the inequality ab < a log a + e b , 
a, b > 0, which leads to the variational formula |26{] 

#K.n|<)= su P f(^ N - 1 ^ / >,;l (4-43) 

feC(x,R) 1 ' } 

for the n-site marginals, L > n, N G N 2 . Note that along with the product measure z/ M 
also the canonical measures n L N are permutation invariant, so that we do not have to 
specify the lattice sites in the n-site marginals. Following (13.41) . 

lim - log Z L ^ L = p(fj,) - ii - p- lim h L ^ L {p) {AAA) 

L— >oo i/ L— >oo 

for all pe D hl . Inserting p = M(p) leads to (l3~8l) . □ 



4.4 Remaining proofs of Section 3.3 

The proof of Theorem 13 . 61 uses large deviation results on the asymptotic distribution of 
T^{rf) as L —>■ oo, which we summarize in the following lemma for our purpose. 

Lemma 4.1 Let u±, u>2, ■ ■ ■ G Z be i.i.d. random variables with mean p c and probabil- 
ity mass function v{k) ~ k~ b with b > 2. Then for every p > p c and some constant c, 
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as L — ► oo, 



{jE^ - > (P-Pc)^}) , (4.45) 

= > cLz/ ({^i = [(P-Pc)^]}) • (4-46) 



Proof. The first statement is shown in [|16ll Chapter 1, Corollary 1.1.1 to 1.1.3, the 



second in 112011 . Theorem 2.2. 



Proof of Theorem 13.61 The proof follows closely the one given in Il20h . Theorem 
2.2 and we only sketch the most important steps. Take p e A{ and fx = M(p). To 
establish (13.261) . using Theorem [3TTI it remains to show that for all e > 

YimJ L ^^{M i L >{ Pi -R c M-e)L})=l, {AM) 

where we use the shorthand M l L (r}) = max x£ A L rji{x). According to (|3.22l) the condi- 
tional measure of the inverse event is 

i>®({Mi<C e L,Vi = [ Pi L] 



p^UMl < C e L\) = ^tt-, — • (4.48) 



where C e := pi — R Cji (p) — e. To prove the theorem we show that this expression 
vanishes for L — > oo. We split the event 

{M l L < C € L} = {[L a ] <M l L < C e L} U {M l L < [L a ]} (4.49) 



for some a G (0, 1) which is chosen below. A basic estimate in Q20|] shows that the 
probability of the second event vanishes for all a £ (0, 1), whereas the first one is 
the crucial part. Analogous to the proof of Theorem 13 . 1 L we single out one side that 
contains the maximum amount of particles. The maximum could be on any of the L 
sites, so 



v^)({[L°]<Mi<C e L,Z j L = [p i L}} 



= CL ^ C0 (*)^ 1,W ({^-i = M-*}) ( 4 - 5 °) 

[L°\<k<C e L 

where the constant CgR accounts for over-counting some configurations where sev- 
eral sites have k or more particles. It is derived in detail in [|2o|l. Using monotonicity 
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of v,r we obtain 

< Ml < C e L, SI = [p,L]}) < 

< CLv^{[L°]) ^^({SU > (RM +e)(L - 1)}) . (4.51) 

With (14.451) the right-hand side is of order 

CLv l ^([U]) ujM({n > (RcAp) + e ) L }) = 0(L 3 - b ^) . (4.52) 

On the other hand, the denominator of (14.481) is at least of order L l ~ h due to (14.451) . 
Thus if we choose a E (2/6, 1), (|4.48l) vanishes for L — > oo which finishes the proof. □ 

5 Connection to zero-range processes 

The ensembles studied above arise naturally as stationary measures of zero-range pro- 
cesses showing a condensation transition which has recently attained much interest. 
Condensation transitions in zero-range processes with a single particle species have 
been studied in two cases: For site dependent jump rates of the particles ll2~7l. [28ll the 
condensate is located at the slowest site. This case is closely related to Bose-Einstein 
condensation into the lowest energy level, which has been studied rigorously using 
large deviation techniques (see Ii29ll and references therein). In this paper we consider 
the case of space homogeneous jump rates that induce an effective attraction between 
the particles [4]. Such models have a number of direct applications, such as network 
dynamics or surface growth, and are particularly important in the study of phase sepa- 
ration in related exclusion models (see [Q] and references therein). 

5.1 Definition 

The dynamics of a homogeneous zero-range process with two particle species on a 
finite, periodic lattice A L is defined by the generator 

2 

= EE 9iM*))pi(y-*)(f{n^)-fto)) ■ (5-D 

i=l x,yeA L 

Here gi{r){x)) E [0, oo) is the rate at which site x loses a particle of species i. It 
jumps to site y according to an irreducible probability distribution and the resulting 
configuration is denoted by rf ,x ~* y . We impose <7»(k) = 4=> hi — and thus the 
process is irreducible on each Xl ; n with fixed particle numbers. For finite lattices 
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the generator is defined for all / : X L — ► R, whereas on infinite lattices there are 
restrictions on the state space and the test functions / [0]. 

It has been shown [|8|, |30J] that for every positive weight w : N 2 — > (0, oo) the 



zero-range process with rates 

w(h - l,fc 2 ) , . w(k 1 ,k 2 - 1) 
w(ki,k 2 ) w{ki,k 2 ) 

has stationary product weight i/; L as defined in (12.21) . independent of pj. So indepen- 
dent of reversibility of the process, the canonical (|2.4I) and grand-canonical measures 
(12.61) are stationary and the results of Section 3 apply to the long-time behaviour of 
such processes. Thus our analysis on the static phase diagram, adopted from equilib- 
rium statistical mechanics, applies also to non-equilibrium zero-range processes. On 
the other hand, dynamic quantities such as dynamic critical exponents or two-time 
correlation functions certainly depend on reversibility. 

Note that (15.21) induces a relation between the rates gx and go, and not every two- 
species zero-range process has stationary product measures [8i |30J]. This is in con- 



trast to single species systems, which always have product measures with w(k) = 
Yli=i m m i s case, dD^ = {fi c } and dD p = {p c } consist only of single points, 

and condensation has been directly related to the asymptotic behaviour of the jump 
rate fl]. If g(k) decays slower than a + 2/k as k — » oo for some a > 0, then p c < oo. 

5.2 Generic examples 

In single- species zero-range processes, the asymptotic decay of the jump rate induces 
an effective attraction between the particles and results in convergence of R(fi) on 
dD^ The same idea was used in (si] to study an example of a two-species zero-range 



process with rates 

n \ nn \ fh{ki + 2)\ k2 / b\ 



3 (k)=0(A; a )(l + ^) , (5.3) 

where 0(0) = and 9(k) = 1 for k > 1. The generic feature is that the rate of species 
2 particles depends only on the presence of species 1 particles, and g% is then chosen 
to fulfill (|5.2I) . This corresponds to the stationary weight (13.201) 

k\\ (k\ + l\ fc 2 
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which we already used as an example in Section 3. For ji 2 < the grand-canonical 
partition function contains a geometric series and can be partially summed, 



fcl! ^ /Jfei + l\*a 



fci=0 v JKl fc 2 =o 1 



yV*. l±h hL_ (55) 



fci=0 



So domz and thus int is a rectangle with int = {/x | /xi, /x 2 < 0}. The parts of 
the boundary that belong to depend on the parameter b, resulting in different phase 
diagrams, as illustrated in Figures [Hand |2| in Section 3. For /ii = the factor 

7z = fcil/ TT (i + b) ~ fcr 6 as fci -> oo , (5.6) 

(1 + % / fi 

in (15.51) determines the convergence properties of R(/x) on 

fe>3 = {n | /xi,/i 2 < 0} , 

3>fe>2 D m = {ax|//i,0 2 <O}\{(O,O)}, 

6<2 L>„ = {^j/ix < 0,/ia < 0} . (5.7) 

The boundary between the phase regions D p and A 2 \ Ai for /i 2 = can be cal- 
culated explicitly as #2(^1, 0) = 1 + -Ri(/ii, 0), /xi < 0, whereas the other boundary 
for /ii = is only given implicitly by i?;(0,/i 2 ) = 9^(0, /i 2 ) (see [8] for more de- 
tails). On top of the stationary phase diagram discussed here, the relaxation dynamics 
of this zero-range process shows an interesting coarsening phenomenon, which has 
been analyzed in ll23ll . 

In the following we consider two other examples which have not been studied 
before. The first one, demonstrating that does not have to be a rectangle, is 

<7i(k) = 9(h) (t^) 2 (1 + Wh) , 92(k) = Y ^-. (5.8) 

Here g 2 oc fc 2 so particles of the second species move independently but are slowed 
down by the presence of species 1 particles, and again g\ is chosen to fulfill (15.21) . The 
corresponding stationary weight is 

w(k) - hl {h + 1)fc2 (5 9) 

W[k) -(l + b) kl fc 2 ! • (5 - 9) 
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Also in this model the partition function can be partially summed and written as a 
hypergeometric function 
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Z ( M ) = e eM2 y - kl \ e^+^Ofci = e ^\F x (l, 1; 1 + 6; e^+^J . (5.10) 

So _D M = |p|e M2 < 0} is a closed set for b > 2 due to (15.61) . Since d-D^ is curved, 
the only non-empty condensing phase region is A\ D A 2 . The resulting phase diagram 
is shown in Figure [3l As in Figure [H the dashed lines point in the normal directions 
of dD^ and in the density plane they determine the background density R c (p). In 
the above example with rates (15.31) these lines are actually uniquely determined by the 
phase boundaries alone (see Figured! right), whereas for this new example they have 
to be fixed via the normal vectors || (1, e^ 2 ) of dD^. Thus Theorem 13.31 implies 
that the background densities fulfill 

ei^M !otMp€AinA2 . (5 . n) 

Pi-Rc,i{p) 1 

Using d5~9l) and (15.101) it is easy to see that R 2 (p.) = (l + #i (/■*)) C 2 for all p G L> M . 
Furthermore 

using standard expansions of the hypergeometric function 2 F\ which are summarized 
e.g. in pO. Taken together, this implies that for p E Aid A 2 

By coincidence, the lines {p | R c (p) = (j^, Pc,2)} converge in the point (—1, 0) for 
all p c ,2 > 0, as is shown in Figure [3] on the right. 

If w is composed of several parts, one can produce various kinds of phase diagrams. 
For example if we add w(k 2 ,ki) to (15.91 ) we get the symmetrized version 

,., fa! (fa + fa! (fa + ljfa 

m(k) ^ (l +% , fa! + (l + % , fa! • (5 ' 14) 

The domain is then given by the intersection of from (|5.9I) with its symmetric 
counterpart, i.e. 

D^= {p\e^ 2 +fn,e^ +fi 2 <0} . (5.15) 

This is illustrated in Figure |4] together with the phase diagram, where phase region 
A\ fl A 2 now shows two different kinds of behaviour of the function R c (p). 



29 



0.2 


-0.2 
/i 2 -0.4 
-0.6 
-0.8 



■1.4-1.2 -1 



-0.8 -0.6 -0.4 -0.2 
Ml 



4 






3 






P2 2 






1 




A in A 2 



-1 



Pi 



(b-2)- 1 



Figure 3: and phase diagram for the zero-range process with rates (15.81 ) for 6 = 4. Dashed 
lines on the left denote the normal directions to dD^ and on the right they determine the 
background density R c (p). 
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In this way one can find zero-range processes exhibiting all kinds of phase dia- 
grams. However, these models are often artificial since the jump rates are very com- 
plicated due to the constraint (15.21) . in particular for the last example. On the other 
hand, simple rates may lead to zero-range processes for which the stationary distribu- 
tion is unknown and not of product form. For such models, a recent study revealed the 



possibility of a discontinuous condensation transition 113 ill . 



5.3 Further remarks 

In section 4.3 we noted that the regularity condition (|3.5I) on the exponential tail of 
w can be relaxed considerably. However, for the results of Section 3 to hold one has 
to assume some regularity of w. Consider for example a single- species zero-range 
process with stationary weight 

»<*) = {£;j^. < 5 -'« 

for some set / = {ii, i 2 , . . .} C N. If |/| = oo then = (— oo, 0]. The proof of 
Theorem 13 .31 works as long as sup JgN (i J+ i — ij) < oo. This is not a purely technical 
condition, because if it is violated, e.g. for ij = 2 J , one does not expect the condensate 
to be stable, since it cannot fluctuate in size. Such a stationary weight leads to jump 
rates with exponentially growing variation as k — > oo, so Monte-Carlo simulations for 
such processes are not feasible. But the behaviour of such irregular processes is in 
general only of limited interest. 

Following the results in Section 3.3, in phase region At PI A 2 we expect the conden- 
sate of each species to concentrate on a single lattice site. Moreover, in the examples in 
Section 5.2 these two condensates are expected to be on the same lattice site, since g\ 
is a decreasing function of k 2 and vice versa, inducing an effective attraction between 



the condensates. Indeed this is what is found in simulations [23]. If both species are 
independent the stationary weight factorizes, i.e. w(k) = Wi(ki) 102(^2), the conden- 
sates do not interact and both have independent random positions. On the other hand, 
if g 1 is increasing in k 2 and vice versa, the condensates repel each other and are not 
found on the same site. In general, whenever the species are coupled, the presence 
of a condensate of one species influences the distribution of the other species on that 
site, also if only one species condenses. This effect is important for the analysis of the 



coarsening behaviour for two species systems and is studied heuristically in [23]. 

All results of Section 3 only address the stationary distribution of zero-range pro- 
cesses. Apart from studying the coarsening dynamics, the ergodic behaviour of the 
system on an infinite lattice A is a dynamical question which is expected to be closely 
related to the stationary results. Starting with a homogeneous distribution /x p (0) with 
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density p = (77(2;))^ , Q >, £ G A at time £ = 0, we expect as t — > 00 

M^MOje^^,, i.e. (/} Mp(t) - (/}^ (p) (5.17) 

for bounded cylinder test functions / G C&(X, M). Although M (see Proposition 
12.31) is the same function as for the stationary results, (15.171) is a statement about the 
dynamics of a zero-range process (15.11) and requires a completely different analysis. 
Such ergodic results exist for attractive single species systems with non-decreasing 
jump rates g(k), where one can use coupling techniques. These are not applicable in 
case of condensation, since then even the single- species process is not attractive, and 
we are not aware of other results in this direction. 



6 Summary 

In this paper we adapt the theory of the equivalence of ensembles to study phase sepa- 
ration in particle systems with unbounded local state space. Our results cover conden- 
sation transitions in (multi-species) zero-range processes which are currently of partic- 
ular interest, and are the main motivation for this study. We use the method of specific 
relative entropy, previously applied to systems with bounded Hamiltonians, which in 
our case involves large deviations and multivariate local limit theorems of subexpo- 
nential distributions. We derive the phase diagram for the condensation and explain its 
connection to the mode of convergence in the equivalence of ensembles, generalizing 
previous results for the non-condensing case. Condensation is shown to be a contin- 
uous phase transition, where the mechanism is different from systems with bounded 
Hamiltonian, and can be characterized by convergence properties of the Gibbs free en- 
ergy on the boundary of its domain of definition. The analysis also involves interesting 
properties of the boundary behaviour of multivariate power series. 

For simplicity of presentation we formulate our results not in the most general set- 
ting, but focus on a particular case which captures the basic novelties of the paper with 
respect to previous work, and is closely related to the main application to zero-range 
processes. A generalization to any number of particle species, each having arbitrary 
discrete state space is straightforward, as long as the stationary measures are of product 
form. Since the method of specific relative entropy only makes use of permutation in- 
variance, an extension to non-product measures along the lines of [|12ll is possible, but 
requires a substantial amount of work. Single- species processes leading to measures 
with nearest-neighbour Hamiltonians have recently been investigated non-rigorously 



fl32H . The result for the equivalence of ensembles is expected to be the same, but the 
structure of the condensate should be different due to spatial correlations. Another 
open point is a rigorous result on the structure of the condensate for more than one 
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species, which involves large deviations for multivariate subexponential distributions. 
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A Results from convex analysis 

In the following we summarize a few results from convex analysis which are used in 
the paper, taken from [|33|l . 

For a function / :K"->1=MU {— oo, 00} we denote by 

dom / = {x G M" I \f(x)\ < 00} (A.l) 

the domain of /. / is called proper if dom / 7^ and f(x) > —00 for all x G M. n . 
For example any function / : D — > R defined on some D C K" can be extended 
to a proper function by setting / = 00 on IR n \ D. To simplify some points we will 
often concentrate on proper functions / : W 1 — > R U {00} in the following. A proper 
function / is convex if for all x, y G R n and r G (0, 1) 

/((I _ T ) x + Ty ) < (1 _ T )/( x ) + Tf ( y ) . (A.2) 

Note that for x G" dom / or y G" dom / (IA.2I) holds trivially, and in particular it 
implies that dom / has to be a convex set. / is strictly convex if (IA.2I) holds with strict 
inequality. We denote by 

argmin/ = Ix G dom / I f(x) = inf f(x)} (A.3) 

iGdom / 

the set of minimizers of /. 

Theorem A.l Let f : IR n — > R be convex. Then argmin/ is a convex subset o/R n . If 
f has a local minimum in x G dom / then x G argmin/, i. e. every local minimum is 
a global minimum. If f is strictly convex argmin/ is either a singleton or empty. 
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Proof. See (13311 . Theorem 2.6. 



/ : R n — > R is called lower semicontinuous at x if 

liminf /(?/):= lim ( inf f(y)) = f(x), (A.4) 

and /ower semicontinuous if this holds for every x G R n . We denote by 

l ev < Q / : = {i e f(x) < a} and lev =Q / := {x G R n \ f(x) = a} (A.5) 
level sets and level lines of a proper function /. 

Theorem A.2 / is lower semicontinuous if and only if the level sets lev< Q / are closed 
in R™ for all a el. If f is convex, then the level sets \ev< a f are convex. 

Proof. See [33]], Theorem 1.6. and Proposition 2.7. 

For a proper convex function, the subgradient at x G dom / is given by 

Sf(x) := {v G R n | f(y) > f(x) +v-(y-x) for all y G dom /} . (A.6) 

If / is differentiable in x, then Sf(x) = { V f (x ) }- There is a more general definition 
of subgradients for non-convex functions (see 11330 , Definition 8.3), which we omit 
since we do not make use of it. It is consistent with (1A.6I) for convex functions as 
proved in $3$ \. Proposition 8.12. 

For any function / : R n — ► R, the convex conjugate function /* : R™ — > R is given by 

f*(v) = sup (v ■ x - f(x)) , (A.7) 

x£M. n 

and the mapping / i— > /* is called the Legendre-Fenchel transform. 

Theorem A.3 Let f : R n — > R a proper convex function. Then f* is proper, convex 
and lower semicontinuous. If f is lower semicontinuous in x G dom /, then 

v e 6 f(x) o x G Sf*(v) o f(x)+f*(v)=v-x. (A.8) 



Proof. See B3h . Theorem 11.1 and Proposition 11.3. 



B Construction of ensembles on a common space 

To formulate the convergence result (|3.7I) we need to define the sequences of canonical 
measures and the grand-canonical product measure on a common measurable space. 
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Take A L C A L+1 and A L | A as L — ► oo, i.e. A L C A for all L and for all 
x E A there exists L > 1 such that a; G A L . Set X = N A and let ^4 be the c-algebra 
induced by the product topology on X, which is generated by the set of all cylinder 
configurations. We can identify Xl by the set of cylinder configurations on A^, which 
are 

r)= {CeX\ C{x) = rj(x) for all x G A L } , (B.l) 

so that A = a(Xi,X 2 , . . .). The family directly extends to a product measure z/ M 
on (X, A) such that f M (r?) = ^(f?) for all 77 G Xl, L > 1. On the other hand, 7Tl,-n 
is a measure on (X, <j(Xl)) where <j(X£) is the cr-algebra generated by Xl, i.e. the 
smallest a-algebra that contains all cylinder configurations 77 G Xl- Then Kl,is ls a ^ so 
a measure on (X, <r(X fc )) for all k < L, since A L C A L+1 implies o"(X L ) C a(X L+1 ). 
For a cylinder function / : X — ► M there exists a fixed n G N, such that / is er(X n )- 
measurable, i.e. / depends only on coordinates in A n . Then we take (X, a(X n )) as 
the measurable space on which ttl,n and are defined for L > n. This construction 
is sufficient to make sense of (13.71) . since it only addresses the limit L — > 00. 
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